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The statistical mechanics of the flux-line lattice in extreme type-II superconductors is studied 
within the framework of the uniformly frustrated anisotropic three-dimensional XF-model. It is 
assumed that the externally applied magnetic field is low enough to invalidate the lowest Landau- 
level approach to the problem. A finite-field counterpart of an Onsager vortex-loop transition in 
extreme type-II superconductors renders the vortex liquid phase-incoherent when the Abrikosov 
vortex lattice undergoes a first order melting transition. For the magnetic fields considered in this 
paper, corresponding to filling fractions / given by 1/f = 12, 14, 16, 20, 25, 32, 48, 64, 72, 84, 96, 112, 
and 128, the vortex liquid phase is not describable as a liquid of well-defined field-induced vortex 
lines. This is due to the proliferation of thermally induced closed vortex-loops with diameters of 
order the magnetic length in the problem, resulting in a "percolation transition" driven by non-field 
induced vortices also transverse to the direction of the applied magnetic field. This immediately 
triggers flux-line lattice melting and loss of phase-coherence along the direction of the magnetic field. 
Due to this mechanism, the field induced fiux lines loose their line tension in the liquid phase, and 
cannot be considered to be directed or well defined. In a non-relativistic 2D boson-analogy picture, 
this latter feature would correspond to a vanishing mass of the bosons. Scaling functions for the 
specific heat are calculated in zero and finite magnetic field. From this we conclude that the critical 
region is of order 10% of Tc for a mass-anisotropy M^/M — 3, and increases with increasing 
mass-anisotropy. The entropy jump at the melting transition is calculated in two ways as a function 
of magnetic field for a mass-ansitropy slightly lower than that in YBCO, namely with and without 
a T-dependent prefactor in the Hamiltonian originating at the microscopic level and surfacing in 
coarse grained theories such as the one considered in this paper. In the first case, it is found to be 
AS = O.lfcs per pancake-vortex, roughly independent of the magnetic field for the filling fractions 
considered here. In the second case, we find an enhancement of AS by a factor which is less than 
2, increasing slightly with decreasing magnetic field. This is still lower than experimental values of 
AS ~ 0.4fcs found experimentally for YBCO using calorimetric methods. We attribute this to the 
slightly lower mass-anisotropy used in our simulations. 

Pacs-numbers: 74.20.De, 74.25.Dw, 74.25.Ha,74.60.Ec 



I. INTRODUCTION 

The physics of vortex matter represents a new field 
of research, which has opened up after the discovery 
of large fluctuation effects in the extreme type-II high- 
Tc pSuperconductprs. In particular, the work of Gammel 
al.U and NelsonEI were important milestones in the field, 
suggesting for the first time that the Abrikosov vortex 
lattice might melt well below the zero-field critical tem- 
perature. Tfjiis extension of ideas originally proposed by 
EilenbergerH, has proved to be fruitful. The melting of 
the fiux-line lattice (FLL) well below the upper critical 
field crossover line due to large anisotropy combinedpSvith 
extreme nonlocality of vortex-vortex interactionsH, as 
well as its first order character, are now well established 
both on theoretical and experimental grounds. Nonethe- 
less, the nature of the molten phase remains tantaliz- 
ingly elusive. Obviously, there exists a 2D boson anal- 
ogy picture of the low-temperature phase of the FLL, 
where the corresponding boson system is an insulating 
one. This is nothing but Abrikosov's mean-field solution 
to the problemO. There is, however, mounting analyt- 



ical and numerical evidence that a similar, intuitively 
appealing, picture of the molten phase via a 2D non- 
relativistic supedfrnd boson systema, may need a sub- 
stantial revisions E3. As far as experimental results are 
concerned, the situation also appears quite intriguingtll. 
At issue here is to what extent the well-defined flux lines 
of the low-temperature FLL phase retain their integrity 
in the high-temperature molten phase. In this paper, we 
address the issue of the character of the vortex liquid via 
extensive Monte-Carlo simulations of the uniformly frus- 
trated three-dimensional XY {3DXY) model. 
An issue of fundamental importance is whether or not a 
vortex liquid (the molten phase of the FLL) is a super- 
conductor or not when the vortex system is not pinned. 
It is clear that the superfluid response to a current ap- 
plied transversely to the field induced vortices is zero at 
all temperatures, provided that pinning is absent. In 
the perfect FLL there will however be a superfluid re- 
sponse to a current applied parallel to the magnetic fleld. 
There remains the possibility that a finite superfluid re- 
sponse might remain even in the liquid phase for this 
geometry, and that it may vanish only deep in the vortex 
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liquid as suggested originally by Feigel'man et aZ.El. A 
main point of this paper is to show that for the mod- 
erate mass-anisotropy and large range of magnetic fields 
considered here, this in fact does not appear to be the 
case. Moreover, the reason that this is so has important 
consequences for the physical picture of the vortex liquid 
phase. Below, we give a summary of the main results of 
this paper. 

The specific heat is calculated for filling fractions / — 
1/12, 1/128,0, and the results are found to be in good 
agpeement with the experimental results of Salapipn et 
alx2 and more rtjcent experiments of Roulin et a/.til, and 
Schilling et altB'lB. In particular, we find that the near- 
logarithmic singularity in zero magnetic field (the spe- 
cific heat exponent a — —0.007) is converted to a broad 
crossover with a peak value suppressed rapidly compared 
to the zero field case. This crossover defines, somewhat 
arbitrarily, the upper critical magnetic field. The entropy 
associated with the suppression of the specific heat at the 
crossover is partly compensated by the appearance of a 
(5-function anomaly in the specific heat at temperatures 
well below the zero-field critical temperature. Such a S- 
function peak is identified unambiguously with the first 
order melting transition of the FLL. The field depen- 
dence of the temperature at which this appears defines 
the melting line of the FLL in the B — T-phase diagram 
of the superconductor. The field dependendence of the 
entropy of the transition is found to be AS ^ B, con- 
sistent with the suppression of the broad main crossover 
peak in the specific heat. This implies that the entropy 
per vortex per layer is essentially field independent in the 
field range and anisotropy range investigate d here. We 
emphasize that the anisotropy considered, yjM^jM = 3, 
is moderate and somewhat smaller than what is found 
in YBCO. Nonetheless, the qualitative aspects-ef our re- 
sults conform well with those found in YBC0ll3. 
We find that in the field regime we have considered, i.e. 
fields down to i? ~ 1 — 5T, the melting of the vortex lat- 
tice is triggered by a proliferation of thermally induced 
closed vortex rings of order the magnetic length of the 
system. This immediately leads to a "percolation" of 
vortex loops traversing the entire system in any given 
direction, and in particular in a direction perpendicular 
to the applied magnetic field. Hence, flux lines which 
are field induced, will traverse the entire system as it 
weaves its way from the bottom to the top of the sys- 
tem. In technical terms, in a simulation one needs to 
apply periodic boundary conditions at least once in the 
(x, ?/)-directions before a flux line starting at the center 
of the bottom layer has reached the top layer. Therefore, 
it does not make sense to view the vortex liquid phase as 
a collection of well defined field-induced flux lines. The 
above picture effectively means that the flux-line tension 
has vanished in the liquid phase. Within the 2D boson 
analogy picture, an equivalent statement would be that 
the boson uLass, which is the analog of the line tension, 
has vanishedtJ. 

Scaling functions for the specific heat are calculated, in 



zero field as well as in finite field. From the regime 
where scaling is found in zero field, we find the width 
of the critical region to be jT — Td/Tc w 0.1 for a mass- 
anisotropy ratio r = y/ Mz/M = 3. This is consider- 
ably wider than what one would naively obtain using the 
Ginzburg-criterion. Moreover, this is slightly wider than 
what we have found for the isotropic case. The width of 
the critical region therefore increases slightly with mass 
anisotropy, for the moderate values of F we have con- 
sidered. Vortex loops are expected to be particularly 
important for the statistical mechanics of the FLL, pro- 
vided that the melting line is found in the proximity of 
the critical region. From the obtained melting curve and 
the width of the critical region, we find that the melting 
curve crosses the critical region curve at a field of order 
_B = IT for F = 3. Below this field, vortex loops will 
completely dominate the physics at the melting transi- 
tion. 

This paper is organized as follows. In section II we 
present the model, the approximations involved and the 
physical quantities considered, as well as updating pro- 
cedure and the parameters used in our Monte Carlo sim- 
ulations. In section III we present and discuss detailed 
results for the filling fraction / = 1/20. In section IV 
we present results in a broad range of filling fractions 
1// S [12, .., 128]. Section V presents the conclusions of 
this paper. 



II. THE MODEL 

The phenomenological model considered in this pa- 
per for the high-Tc cuprates, is the uniformly frus- 
trated 3 jiinieiasipnal anisotropic XY {3DXY) model on 
a lattice JlTt^a'El, defined by the Hamiltonian 

Hmr)}) = - E cos[V^^(r)-A^(r)], (1) 

r,fj,—x^y,z 

where 9 is the local phase of the superconducting complex 
order parameter and V is a lattice derivative. Further- 
more, the coupling energy along the /i-axis, J^, is defined 



by 



Jx Jy — 



Here $0 is the flux quantum, S^ab is the superconduct- 
ing coherence length within the CuO-planes, and d is the 
distance between two CuO-layers in adjacent unit cells. 
Furthermore, Xab and Ac are the magnetic penetration 
lengths in the CuO planes, and along the crystals c-axis, 
respectively. In Eq. |^, is related to the quenched 
vector potential A^p by 
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where is the unit vector along the fi-axis. This 3DXY 
model is dual to the i-apisotropic London model in the 
limit of (Aab, Ac) — > ocEil. This limit should be taken with 
the understanding that the coupling energies J± and Jz 
are maintained finite. When (Aa&, Ac) — > oo, gauge fluc- 
tuations are completely suppressed, leaving a quenched 
vector potential and uniform magnetic induction. Thus, 
the 3DXY model should give an adequate description 
of the physics of extreme type-II single crystal supercon- 
ductors in the field regime Bd << B « Bc2- The 
condition Bd « B ensures that the contribution to 
the magnetic induction from individual flux lines over- 
lap strongly giving a uniform magnetic induction. The 
condition B « Bc2 ensures that details of the internal 
structure of the vortex cores are not essential. Moreover, 
the 3DXY model should give an adequate description of 
the physics of extreme type-II single crystal supercon- 
ductors in zero-magnetic field when gauge fluctuations 
are not important. 

In this paper we consider simple tetragonal systems 
with dimensions = Ly = L± and L^- The co- 
ordinate (x, y, z)-axes are taken to be parallel to the 
crystal (a,b,c)-axes, respectively. We measure the in- 
plane length scales {x,y, L^, Ly) in units of ^ab and the 
length scales along the z-axis (z, L^) i units of d. Our 
unit cell is a simple tetragonal system with dimensions 
S'x — Gy — S^ab, Cz — d. Periodic boundary conditions are 
used in all directions throughout. 



A. The internal energy and specific heat 

The specific heat per site C is obtained using the stan- 
dard fiuctuation formula 



1 <H'> -<H>' 



[kBTf 



(2) 



where ks is Boltzmann's constant. Results for most tem- 
peratures are checked for consistency by differentiating 
the results for the internal energy with respect to tem- 
perature. To estimate the latent heat (entropy jump) at 
a first order phase transition, we consider the internal 
energy per site E, 



E = 



1 



L I L, 



<H> 



(3) 



This expression holds as long as we do not include any 
T-dependence in the Hamiltonian. Such a T-dependence 
could conceivably arise in effective coarse-grained theo.- 
ries such as the GL-theory, as first pointed out in Ref.cS. 
At a first order phase transition there is a discontinuity 
in the internal energy per site AE associated with coexis- 
tence of the Abrikosov FLL phase and the vortex liquid. 
ThisJii turn gives rise to a J-function peak in the specific 
heatErLl. For the melting transition of a vortex line lat- 
tice, AE is related to the entropy jump per vortex lines 
per layer AS by 
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AE 
fksTn 
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where T„i is the melting temperature and / is the vortex 
lines density defined below in Eq. ^j. For consistency, one 
may also check this result by extracting the entropy jump 
at the melting transition from the scaling ofpthe-height 
of the (5-function anomaly in the specific heatc3i3'Ll 



C-const+-{^\ . (5) 



B. Scaling functions for the specific heat 

An issue of principle importance is whether or not 
closed thermally induced vortex loops, i.e. the criti- 
cal fluctuations, will influence the melting of the FLL. 
Naively, one expects this to be the case provided that 
the melting temperature Tm{B) is within the critical re- 
gion of the zero-fleld transition. It is therefore a matter 
of interest to establish the width of the critical region 
of the anisotropic 3DXY-niode\. To this end, we con- 
sider flnite-size scaling of the speciflc heat. The region 
of data-collapse of the speciflc heat evaluated at various 
system sizes identifles the width of the critical region. 
It is also of interest to find the extension of this critical 
region to finite fields, i.e. the width of the crossover re- 
gion around the upper critical field, and how it depends 
on the anisotropy ratio F. If the melting line is located 
within this crossover-ijegion or close to it, one expects a 
vortex-loop "blowout"cli3 to dominate the physics at the 
melting transition, analagous to what-was suggested to 
happen in superfluid He'^ by OnsageiO. This is partic- 
ularly important at very low magnetic flelds_iike those 
considered in the experiments of Zeldov et aZ.Ell. 
For the zero- field finite-size scaling of the specific heat, 
we have used cubic samples L x L x L with L = 
16,32,48,64,72,96 to avoid spurious geometric effects. 
To investigate the possibility that the width of the criti- 
cal region may depend on anisotropy, we have considered 
the two cases F = 1 and F — 3. 

The finite-size scaling function for the specific heat may 
in general be obtained in stanci ard fashion from the sin- 
gular part of the free energy asE£lc3 



C(i,L) = L"/'' <^±{\t\L^^''), 
or equivalently 

G [t, oo) 

where |t| ^ \T - T^\/T^, and where ^±{\t\L^''') and 
G±{\t\L^^") are analytic functions of their arguments. 
Here, a is the specific heat critical exponent, and v 
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is the critical exponent of the superconducting correla- 
tion length, hyperscaling yields a = 2 — Dv in S-i-D- 
dimensional system. As discussed in detail in Ref.E£l, a 
more convenient scaling form for numerical purposes is 
given by 



C{t,L) - C(0,oo) 
C(0,L) - C(0,oo) 



(6) 



We will use this scaling form to determine the width of 
the critical region. 

In the presence of a general field X with scaling dimen- 
sion X ^ where f is the correlation length ^ ~ l^l""; 
we have 



c(i,x) = |t|-" g±[x\t\ 



-\u\ 



Subtracting out the zero-field part, and introducing y 
X \t\-^'' and AC(t, X) = C{t, X) - C{t, 0), we find 



AC(i,X) 



Q±{y)-Q±{Q) 

Q±{y)-Q±{Q) 



n±iy). 



Choosing X = B, the induction, implies that the scaling 
dimension A = 2 when gauge-fiuctuations are suppressed, 
as is the case in the uniformly frustrated 3DXY-m.ode\. 
Under such circumstances, the induction B will not ac- 
quire anomalous scaling. Hence, we obtain 



in three dimensions, we have \t\ ^ B^^^, implying that 
the crossover region around the upper critical on the low- 
temperature side has a positive curvature in the B — T 
phase-diagram, which is also true for the melting curve, 
for which we have \t\i\i ~ B'', with t] ~ 2/3. The widen- 
ing of the crossover region is of course consistent with a 
broadening of the remains of the zero-field anomaly in 
the specific heat, to be calculated below. 
The crossover curve B{T < Tc) has a more rapid increase 
as a function of Tc — T than the melting curve, recall the 
exponents 3/4 and 2/3, respectively. Due to the finite 
width of the zero-field critical region, there should then 
be a field regime where either the melting curve and the 
crossover curve intersect, or where the crossover curve 
is to the left of the melting curve in the B — T phase- 
diagram. This depends on the width of the zero-field crit- 
ical regime. Given the size of this regime, \t\ < 0.1, the 
former scenario appears to us to be the more likely one, 
and this is also what we find in our simulations. Hence 
critical fiuctuations, i.e. thermally induced vortex loops, 
should substantially influence the FLL melting in a finite 
regime of magnetic fields. From our simulations, to be 
presented below, we estimate the relevant field regime to 
be of order — IT in an extreme type-H superconductor 
with r = 3. 



C. The helicity modulus 



^a/2u ^c{t,B) 



y 
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As a probe of global superconducting phase-coherence, 
G±iy) — g±{0) = 'H±{y), (7) consider various helicity moduli T^;, Tj, and T^. The 

helicity modulus along the /i-direction is defined as 
the second derivative of the free energy with respect to a 
global phase twist along the /i-directior£j, explicitly we 
obtain for the anisotropic uniformly frustrated 3DXY- 
model 



where y = B\t\~'^'^ . We will use the above scaling 
form Eq. ^ to determine the width of the crossover 
region around the upper critical field, as a function of 
B. Note, however, that this scaling form is not specific 
to the 3£'A'y-model. By plotting appropriate ratios of 
temperature- and field-derivatives of these scaling func- 
tions, one may extract directly the critical exponents of 
the system as the slopes of the quantitities being plot- 
ted, see for instance the very detailed analysis of this by 
Schilling et alS3. It is conceivable that such a proce- 
dure would yield a curve with a kink in it when i > 0, 
as claimed to be observed by Schilling et al. This in 
itself does not invalidate the 3Z)A'F-scaling of high- Tc 
cuprates. Conceivably, it could be due to a crossover 
from an XY fixed point to another fixed point, possibly 
with an anomalously large value of v « 1.5, based on 
magnetization data. (Note that the specific heat data of 
Schilling et al. in fact show an opposite trend, more con- 
sistent with a crossover to a Gaussian fixed point. This 
is to be expected if amplitude fiuctuations of the order 
parameter were to dominate the phase- fiuctuations) . For 
more details, see the discussion below. 
We note immediately that the above implies that \t\ ~ 
^1/21^ for finite fields, i.e. the width of the crossover re- 
gion widens as B increases. Using the estimate v — 2/3 



T—r ( ■^'^ cos[V^6l(r) - A^{r)]{e^ ■ 

^±-t^z \ 

\r,u—x.y,z 



kBTL±Lz 



^ sin[V^6i(r) - A^{r)]{e^ ■ e^) 



_r,i'—x,y,z 



When is finite, the system can carry a supercur- 
rent along the /z-direction. When vanishes, resis- 
tivity along the /i-direction becomes finite. In systems 
with finite applied field along the z-axis, we expect 



for all temperatures in the continuum limit. 



In this case, any applied current along the xy-plane will 
move the unpinned flux lines and dissipate energy. Dis- 
cretization introduces a potentially singular perturbation 
by introducing an artificial pinning potential, the effect of 
which is more serious in a three-dimensional system than 
in a two-dimensional one. In the latter case, the effect 
of the potential may in principle be eutirely avoided by 
considering low enough filling fractionsEZI, whereas this is 
not possible in three dimensions in the thermodynamic 
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limit. The size Lz of systems must therefore be tailored 
to the filling fraction / in order to avoid spurious pinning 
effects. 

The thus introduced pinning potential will, at a low 
enough temperature, pin the flux lines in their positions, 
and cause (Ta;,Tj,) 7^ up to a depinning temperature 
Td- To ensure that this artificially introduced pinning po- 
tential caused by the numerical lattice does not affect the 
FLL melting transition, we should consider systems with 
Td much lower than all other "critical" temperatures of 
interest. Td is controlled piainly by the filling fraction /; 
we have ^ as / ^ 0EZH3. To adequately mimick the 
continuum limit of interest, low enough filling fractions 
must therefore be considered. 



D. The FLL structure function 

To locate the position of the vortex elements we use 
the following procedure: The counterclockwise line inte- 
gral of the gauge-invariant phase-differences around any 
plaquette of the numerical lattice with surface normal 
along the /z-direction must always satisfy 

^j.(r) = 27rK(r)-^), 

>(r) = V,^^(r) -^,(r). 

Here, Ci is the closed path traced out by the links sur- 
rounding an arbitrary plaquette, and v represents the 
Cartesian components of the current in the directions of 
the links which comprise the closed path Ci. Further- 
more, ji/(r) is the current on the link between site r and 
site r -|- e,/, and n^ir) = 0,±1 represents a vortex seg- 
ment penetrating the plaquette enclosed by the path Ci . 
Here, is the vortex lines density along the /^-direction, 
and is given by 



L\Lz 



(8) 



To probe the structural order of the vortex-system, we 
consider the inplane structura function for vortex seg- 
ments within the same planeB, 



5(ki) = 



PL'\ Lz 



E 



\nz{v±,z)e 



In the FLL phase we expect to see a periodic array 
of sharp Bragg peaks in the k^-plane. In the vortex 
liquid phase we expect to see Bragg rings with radius 
fc_L = I'K I and 47r/a^, characteristic of a liquid. Here, 
is the average distance between neighboring vortex 
lines. 



E. Monte Carlo procedure 

The Monte Carlo updating procedure used in this pa- 
per is the following. The numerical lattice is stepped 
through in a systematic manner. At each site a change 
of the local phase of the superconducting condensate is 
attempted by a random amount S [— 7r,7r). The at- 
tempt is accepted or rejected according to the standard 
Metropolis algorithm. 

If the accepted phase change causes the current on a link 
j^(r) to fall exceed the range j^iir) € [— tt, tt), an amount 
d=27r is added to the current such that j^{r) is brought 
back into the primary interval j^if) S [— 7r,7r). An im- 
portant point is that this operation can only generate a 
closed unit vortex loop around the link where the current 
is changed, thereby conserving the net induction of the 
system. No net vorticity is ever introduced by the proce- 
dure, and the procedure also guarantees that no vortex 
line can start or end within the sample. One sweep refers 
to L^Lz attempts to change the phase angle. 
We fix the height of our systems to Lz = 40 and let L± 
vary from 40 to 128 depending on the ffux-line density 
under consideration. In RefsJjfl, it was noted that for sys- 
tems with moderate anisotropy (F ^ 3), finite size effects 
are rather small when the linear dimension of the system 
and the total number of fiux lines exceed ^ 40. Thus, 
we believe that finite size effects will not affect the con- 
clusions in this paper. Likewise, it was observed by the 
same authors that finite-size effects were negligible when 
Lz was increased beyond Lz = 40 for the anisotropy con- 
sidered here, F = 3. This has motivated our choice of 
Lz = 40. 

In this paper we fix the anisotropy parameter F to 



r = , — - 



^ab£.ab 



= 3. 



in most simulations. Occasionally, comparison is made 
for the isotropic case F = 1. The magnetic field B is ap- 
plied along the crystal c-axis, giving a vortex line density 



fx — fy 



0, 



/. = / 



(9) 



The flux-line densities / considered are: 1// = 12 (48), 
14 (56), 16 (48), 20 (40), 25 (50), 32 (64), 48 (48), 64 (64), 
72 (72), 84 (84), 96 (96), 112 (112), 128 (128), 00 (64). 
The numbers in the parentheses denote L± for the corre- 
sponding vortex line density. Note that we have chosen, 
with our gauge, L± for each filling fraction / in such a 
way that we ensure that an integer number of magnetic 
Brillouin-zones will fit on the reciprocal lattice of each 
system, enabling us to use periodic boundary conditions 
in the a;, y-directions. As will be observed, L± is an in- 
teger number of 1// in each case. 

The value of / is prescribed by loading the following 
phase-difference pattern onto the numerical lattice (using 
Landau gauge) a system with 9{r) = for all r. 
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Ay{x,y,z) = 2TTfx. 

The system is then heated to a temperature well above 
any transition/crossover temperatures of interest, at 
which point slow cooling is started. The filling frac- 
tion / is conserved by our Monte Carlo procedure, and 
the moves are carried out on the gauge-invariant phase- 
differences on each link. Note that we do not need to as- 
sume any ground state configuration by this procedure. 
Extremely long simulations, typically 4 x 10^ — 6 x 10^ 
sweeps, are however required in order to capture the cor- 
rect physics at the FLL melting transition and to reveal 
any (5-function anomalies in the specific heat at the melt- 
ing transition, particularly at low filling fractions. 

III. RESULTS, F = 1/20 

A. FLL melting and phase-coherence 

To identify the possible different phases and phase 
transition(s)/crossover(s) in a system with finite flux-line 
density, we first concentrate on results for the system 
/ = 1/20. Similar results are found in all other finite 
flux-line densities considered in this paper, to be detailed 
below. We have measured temperatures in units where 
kB = 1. 

Fig. |l| shows the specific heat per site C, the helicity 
modulus along the applied field direction T^, the helic- 
ity modulus perpendicular to the applied field direction 
Tx, and the inplane structure function S'(k^), as func- 
tions of temperature. Fig. |l| shows that the inplane struc- 
ture function ^(kj^ — 27r/5,7r/4) has a sharp drop from 
0.2 to precisely at Tm — 0.531Jl, indicating that the 
FLL melts at Tm in a first order phase transition. For 
a more global view. Fig. || illustrates the density plot of 
5(k^) for kx, ky g [— tt, tt] at four different temperatures; 
T/J± = 0.450, 0.530, 0.531, 0.532. It is clearly seen that 
the periodic array of sharp Bragg-peaks is converted into 
a ring precisely at Tm = 0.531J^, within a narrow tem- 
perature region of AT — 0.001 Jj^ around Tm- 

To clarify whether the phase coherence along the di- 
rection of the applied magnetic field is finite in the vortex 
liquid phase, we consider the helicity modulus along z- 
axis, Tz ■ In Fig.|l] it is clearly seen that shows a sharp 
jump from 0.6 to at Tm precisely where the FLL melts. 
This shows that the FLL melts directly into an incoher- 
ent vortex liquid in a first order phase transition. We will 
return to this important point later, since it has impor- 
tant consequences for the physical picture of the vortex 
liquid phase. The above result is in complete agreement 
with the. work of RefEl using the SDXY-model, the work 
of Ref.L3 using the lowest Landau-level approximation, 
and earlier work by us using the 3D anisotropic Villain 
modela. In all these works, it was found that longitudi- 
nal phase-coherence is lost as soon as the vortex lattice 
melts in the thermodynamic limit. We emphasize that 
opposite conclusions were drawn in earlier work by us 



and othera3't3E3. We believe that this discrepancy may 
be due to one or several of the following three factors: 
i) In earlier work, the system size in the z-direction may 
not have been large enough, particularly for the isotropic 
case, ii) the simulations were not run for a long enough 
time, and iii) the results were obtained upon heating only. 
Our more recent results in Ref.u and in the present paper 
are obtained upon heating and cooling. 
A first order phase transition is manifest in the form 
of a ^- function anomaly in the specific heat. From the 
height of this anomaly one may deduce the latent heat 
of the transition. Figj^ shows that the anomaly occurs 
at Tm — 0.531Jj_, precisely where the structure function 
and the helicity modulus vanish. These results are in 
complete agreement with those of Rcfs.B obtained on the 
uniformly frustrated 3iiXy-model for / — 1/25, as well 
as those found in RefEl using the uniformly frustrated 
3-D anisotropic Villain-model for / = 1/32. Below we 
consider filling fractions down to / = 1/128, finding that 
these results still hold. 

The latent heat of the first order FLL melting transi- 
tion at Tm is obtained from the jump in the internal en- 
ergy shown in Fig. ^, using Eq. |4[ Here, the entropy jump 
per vortex line per layer is estimated to be AS' = Clfc^. 
To obtain the spike in the specific heat we must i) find the 
transition temperature must be located very accurately, 
typically to within a part in 10^ and ii) increase the sim- 
ulation length to at least 6000000 sweeps over the lattice 
for each temperature. The extreme length of the simula- 
tions is necessary to allow the system to switch back and 
forth between the ordered phase and disordered phase at 
the phase transition an adequate number of times, typi- 
cally at least ten times. 



B. Breakdown of the 2D boson analogy 

The specific heat has a broad anomaly at Tbc2 — 
1.05J± 3> Tm, indicating a crossover. This broad 
crossover was unambiguously identified in our previous 
work as the xemains of a zero field Onsager vortex loop 
"blow ouf'tO'El, destroying superconductivity on all length 
scales. However, since the remains of the zero field vor- 
tex loop "blow out" takes place first at Tbc2 ^ Tm, su- 
perconductivity still exists locally in finite domains in 
the incoherent voptex liquid phase, giving strong diamag- 
netic fluctuation£3. Since the global phase coherence in 
all directions is destroyed in the incoherent vortex liq- 
uid phase, the superfluid stiffness is zero in all directions 
in this phase, and any applied current through the sys- 
tem will dissipate energy. Thus, in the incoherent vortex 
liquid phase the system has both finite resistivity in all 
directions, as well as strong diamagnetic fluctuations. 
The numerical lattice is a singular perturbation in a three 
dimensional system, and one may ask whether the flrst 
order FLL melting transition at Tm is affected by the 
artificially introduced pinning potential. To address this 
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issue we consider the helicity modulus along the x-axis, 
Tx- In Fig. |l], it is seen that the helicity modulus along 
x-axisTa; drops to zero already at = 0.1J± <^ T„j. 
From this we conclude tkat above Td the system exhibits 
a "floating solid" phaseED. Thus the FLL melting transi- 
tion at T„i is not affected by the pinning potential caused 
by the numerical lattice. 

Snapshots of the FLL, the incoherent vortex liquid phase, 
and the normal metal phase are shown in Fig. of 
the system / = 1/20 for four temperatures T/J± = 
0.26, 0.50, 0.54, 0.70. For clarity only a part of the 
system; x,y £ [0 : 20],z G [0 : 40], is shown. For the 
system / — 1/20, we have found Tm = 0.53J_l and 
Tbc2 = 1.05J_L. For T = 0.26J_l < T,„, the flux lines 
form a hexagonal lattice. Although there are many ther- 
mally induced defects attached to each flux line, they re- 
main well defined entities. For T = 0.50J^ ^ T,„, though 
the flux lines now fluctuate substantially, they nonethe- 
less remain intact. So does the FLL, as evidenced by the 
results for the structure function, see Figs. |^ and |^. For 
a slightly more elevated temperature T = 0.54Jj_ ~ T^, 
the FLL has melted. A key observation is that, immedi- 
ately upon melting, the flux lines are no longer well de- 
fined entities, there are many intersections between flux 
lines, vortex-loops have proliferated, and there exists at 
least one way to percolate from one side of the sample to 
the opposite side in any direction. Thus, for any given 
direction there always exist at least one "infinite" long 
vortex lines perpendicular to it, and any applied current 
will move these "perpendicular" vortex lines and dissi- 
pate energy. Note that in this picture vortex lines in the 
incoherent vortex liquid cannot bp, described as world 
lines of 2D non-relativistic boson^. Thus, one vortex 
line in the center of the system will meander all the way 
to the boundary surface (with surface normal perpen- 
dicular the applied field direction) and back as a field 
induced flux line weaves its way from the bottom to the 
top of the system. This corresponds to zero flux-line ten- 
sion, and a wandering exponent C of the flux line which 
is C > 1 or, equivalently, zero bosonic mass in the 2D 
boson analogy. 

Note that the 2D quantum boson sjstem we have in mind 
when referring to the work of Ref.l3 is a non-relativistic 
system. The picture we have in mind for the liquid phase 
is more akin to a relativistic 2D quantum boson system, 
where the proliferation of vortex loops and overhangs in 
the flux-lines correspond to vacuum-fluctuations in the 
boson sjjstem. This connection has been nicely exposed 
in Ref.E3. 



IV. RESULTS, 1/F G [12, .., 128] 

A. Structure function S'(kx) 

We show in Fig. |^ the inplane structure function S'(kj^) 
as a function of temperature for several vortex line den- 



sities /; 

1// = 12(kj^ = 57r/12,57r/12), 16(k_L = Stt/S, -7r/3), 
20(kj_ = 27r/5,7r/4), 25(k^ = 67r/25, 97r/25), 32(k^ = 
57r/16,-77r/32), 48(k^ = 7r/6,7r/4), 72(k^ = 
7r/9,-27r/9), 96(k^ = 37r/16, 57r/58), 128(k_L = 
77r/64, 57r/32). For a given vortex line density, 5'(k^) 
with the corresponding value of k^ shows a sharp drop 
from ~ 0.2 to zero defining a field dependent FLL melt- 
ing temperature Tm{f)- This clearly shows that, for all 
values of / considered here, the FLL melts in a first order 
phase transition. For decreasing f, the transition temper- 
ature T„i{f) increase towards Tc as expected, see Fig. |^. 



B. Helicity modulus along the field direction, Tz 

In Fig. H the helicity modulus along the direction of 
the magnetic field is shown as a function of temperature 
for the same set of flux-line densities as for the case of 
the inplane structure function. As / is varied, shows 
a sharp drop towards zero precisely at the corresponding 
FLL melting temperature Tm{f)- Thus, we may conclude 
that for all filling fractions considered, the FLL melts di- 
rectly into an incoherent vortex liquid. The temperature 
region where the vortex liquid and the phase coherence 
along the applied-field coexist, found previously by sev- 
eral authorsolH3'LrEj, is not found for any flux-line density 
/ considered in this paper. We believe that the temper- 
ature regime where the vortex liquid exist with phase 
coherence along the field direction pertains to thin film 
geometries, or are otherwise an artifact too short sim- 
ulations and hysteretic behavior in the heating/cooling 
sequence of the vortex system. The phase-coherent vor- 
tex liquid does not exist in the thermodynamic limit of 
an equilibrium system, at least in systems with moderate 
anisotropy and moderate magnetic induction. The pos- 
sibility of the existence of a very small magnetic field in- 
duction Blower (dependent on the anisotropy -B/ou)er(r)) 
below which the phase coherence along the field direction 
can exist in the vortex liquid is not completely ruled out 
by this work. 

Note that Tm{f) is correlated with the temperature 
where the corresponding start to fall sharply towards 
zero, not the lowest temperature where vanish. One 
may question whether it is correct to take the tempera- 
ture where T z starts to show a sharp drop as the tem- 
perature where phase coherence along the applied field 
vanish. For moderate vortex line densities this poses no 
problem, since the drop in is very sharp. However, 
for / < 1/48, the transition extend over a small tem- 
perature region. By experience, we know that when the 
system size and the number of vortex lines in the system 
increases, the drop in sharpens and the tail in dis- 
appear. We believe therefore that this tail is only a finite 
size effect. 

The main conclusion of the above discussion is that in 
the thermodynamic limit, no phase-coherence exists in 
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the vortex hquid phase. This conclusion is cons;,s, 
with mounting evidence from numerical simulationf 
obtained however only in a limited filling range 1// G 
[25 — 36]. Our present results extend these conclusion to 
much lower filling fractions. 



C. Specific Heat 

For all filling fractions down to / ~ 1/32, we have 
found a (5- function anomaly in the specific heat at T^if), 
indicating a first order phase transition, see Fig |^. For 
smaller filling fractions / < f/48 we find no clear ev- 
idence of a spike in the specific heat. Note that in 
passing from the system with / = 1/32 to the system 
with / — 1/48, the number of field induced flux lines is 
reduced from 128 to 48. We believe that the observed 
"none-existence" of the (5-function anomalies at the FLL 
melting temperature in system with very low flux-line 
densities is attributable to two factors: 1) for system 
with / < 1 /48 we have to few field induced vortex lines 
in our systems, and 2) the contribution to the specific 
heat from the field induced flux lines for these flUing 
fractions is too small compared to the "spin-wave" and 
vortex-loop contributions, to be detected by our simula- 
tions. 

In Fig. we show the specific heat as a function 
of temperature for the same set of corresponding sys- 
tem sizes and flux-line densities as previously used in 
calculating the inplane structure function and the helic- 
ity modulus along the direction of the applied magnetic 
field. For decreasing / the crossover temperature Tbc2(/) 
increases and moves towards the zero-field critical tem- 
perature Tbc2(/ = 0) = Tc. The broad anomaly in 
the specific heat sharpens and the maximum height of 
the cusp increases, evolving smoothly towards the zero- 
field specific heat singularity at Tc- Tbc2 denotes the 
crossover temperature at which the remains of the zero 
field vortex-loop "blowout" takes place. In a finite mag- 
netic field the vortex -loop "blowout" at TBc2{f) causes 
only a crossover and the actual phase transition takes 
place at a lower temperature, Tm{f), where the FLL 
undergoes a first order melting transition triggered by a 
proliferation of vortex loops with diameters at least of 
order the magnetic length in the problem. 

Scaling functions for the specific heat both for zero 
field and finite field are shown in Figs. || and |l^. 
For the zero-field case, using Eq. ^, it is seen from Fig. 
^ that data-collapse is obtained over a wide region out 
to values of the scaling variable \t\V-/'^ > 10 for L > 32 
on the low-temperature side of Tc. Note also that the 
width of the scaling regime is slightly larger for F = 3 
than for r = 1. We expect this trend to persist with 
increasing F; in the extremely case where the layers may 
be considered completely decoupled, i.e. F oo, the 



entire low-temperature regime is known to be criticalEi. 
For smaller L, it appears from our simulations rtlaat we do 
not obtain scaling. Using the value v = 0.66923, we find 
that the width of the critical region is given by \t\ « 0.1. 
The critical scaling of the specific heat is also consid- 
erably better above Tc than below. This is due to the 
fact that vortex loops, i.e. the critical fluctuations, to 
a much larger extent dominate the free energy above Tc 
compared to below Tc. Below Tc there is a non-singular 
contribution to the free energy, and hence specific heat, 
due to spin-wave fluctuations of the local phase of the 
order parameter. 



The scaling function of the specific heat in a finite field, 
given above in Eq. 0, is also calculated for filling frac- 
tions / given by 1// = 12,16,20,25,32,48,72,96,128 
with correponding system sizes identical to those used for 
the structure function above. The anisotropy is F = 3. 
The scaled results are good agceement with the worka^ 
Salamon et alx3, Roulip et a/.llil, and Schilling et alx3. 
Note that while Ref.ll3 makes the point that 3DXY- 
scaling does not appear to describe well the experimental 
results above Tc in a single-crystal YBa2Cu307 in the 
field-range 0.75 — 7T, it does appear to work well below 
Tc- The reason for this may be that for an optimally 
doped compound, the temperature at which a pseudo- 
gap opens up may not be much higher than Tc at which 
phase-coherence is established. The results are therefore 
likely to be influenced by amplitude fluctuations above 
Tc- This should not be the case below Tc. Hence, we be- 
lieve that the field range considered is not the only issue, 
but also that one is observing a crossover from an XY- 
critical point to a Gaussian critical point when increasing 
the temperature above Tc, approaching a mean-field like 
temperature Tmf where preformed pairs start to be dis- 
sociate. An obscuring factor is that the specific heat data 
and magnetization data of Schilling et al. show opposite 
trends in their deviation from XF-scaling. Note that the 
analysis of RefJlj is not specific to the XY-model, the 
scaling forms that are used are quite general. Were the 
temperature scales Tc and Tm f to be well separated, XY 
critical scaling would presumably persist above Tc. This 
would for instance be the case in underdoped cuprateso. 
At any rate, it is the width of the critical region below 
Tc which is of interest in establishing the importance of 
interplay between vortex loops and FLL melting. The 
width of the critical region should increase with under- 
doping, and hence the interplay between vortex loops and 
FLL melting is expected to be morejoronounced when the 
cuprates become more underdopedE3. 
We note that the scaling is better above Tc than below, 
again because non-singular contributions to the free en- 
ergy, in this case also arising from the FLL, contribute 
significantly. The spikes in the finite-field scaling func- 
tion are due to the specific heat anomalies at the FLL 
melting transition. 
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D. Entropy discontinuity at the FLL melting 
transition 

The latent heat, or equivalently the discontinuity in 
entropy at the FLL melting trajpi tian ,, ha-s , been much 



focused on in recent experiments 



InFig.p, 

the entropy discontinuity at the first order FLL melting 
transition is shown as a function of the flux-line den- 
sity. The results obtained using the Hamiltonian in Eq. 
|l| is shown in filled circles. We find that the entropy 
discontinuity per flux line per layer AS{f) O.lks- 
The fact that AS is essentially independent of the ap- 
plied magnetic field, for the moderate anisotropy F = 3 
considered in this paper, is consistent with, the experi- 
mental results obtained by Schilling et al.lla and Roulin 
et alJlil. They found AS{B) - O.S/fcs, independent of 
B. The values of AS ~ O.lks are similar to the values 
found by Hu et al.El. We attribute the difference between 
our values for AS{f) ~ O.lkB and the experimental 
value AS{B) ^ O.Sfcs to the difference in the anisotropy. 
YBCO has an anisotropy F ~ 7, while the anisotropy in 
this paper is F = 3. As shown in our previous paperQ, 
and also by Hu et al.l, the entropy jump at the FLL 
melting transition increases with increasing anisotropy. 

To ensure that the artificial pinning potential intro- 
duced by the numerical mesh does not affect the FLL 
melting transition at Tm(/), we must ensure that the 
helicity modulus perpendicular to the applied field van- 
ishes at a temperature Td{f) significantly below Tm{f)- 
Under such circumstances, the low temperature phase 
for Td{f) < T < Tm{f) is characterized by a "floating 
solid phase", mimicking the continuum limit. In Fig. 
the helicity modulus along x-direction is shown as a 
function of temperature for the same set of / used for 
the specific heat, structure function, and T^. Fig. |l^ 
shows that for each flux-line density considered, T^: van- 
ishes at a temperature Td{f) significantly lower than the 
corresponding FLL melting temperature Tm{f)- Thus, 
we have shown that in all systems considered in this pa- 
per, the depinning crossover at Td{f) does not affect the 
FLL melting at Tm{f) > Although we have not 

shown it explicitly here, we have checked that Ty{T) is 
essentially identioal to T^^, as required by symmetry. 

In recent workta, it was pointed out that calculated 
entropy jumps AS at the melting transition of the 
Abrikosov vortex la.ttice could be brought into agreement 
with experimentsEj by introducing temperature depen- 
dent parameters in the theory, reflecting fluctuations at 
a microscopic level surfacing in coarse grained theories. 
The idea pf. using such a procedure was flrst introduced 
by in RefEa within the lowest Landau level approach to 
the same problem, i.e. the high-field limit. This leads to 
an internal energy 



U{T) =<H > -T < 



dH 
'df 



(10) 



< H >= {l/Z)J2 H expi-H/keT), and Z = 
^ eKp(—H/kBT) is the canonical partition function. For 
a derivation of this result, see Appendix A. 

In extreme type-II superconductors, as modelled by 
the 3£'Ary-model or the London model in the A oo- 
limit, the T-dependence described above appears exclu- 
sively as a prefactor in the Hamiltonian, H = Eo{t)Ho, 
where Hq has no T-dependent prefactors, r = T/Tcmf 
with TcMF a mean-field zero-field transition temperature, 
and EqIt) = [A(0)/A(r)]2. Hq is to be identified with the 
Hamiltonian used in this paper so far. For instance, in 
the two-fiuid model Eq{t) = 1 — r^, while the simplest 
mean- field approximation yields Eq{t) — 1 — t. Using 
the above, we find the internal energy given by 

U=[Eo{r)-T^^] C/o(T'), 



r 



T 



Eoir)- 



(11) 



This leads to an entropy jump at the first-order melting 
transition of the Abrikosov vortex lattice 



AS- 



AU _ 
T 
1 

Eo{t) 



[Eoir) - 
Eoir) 



T 



T 



dEoir) AUojT') 

dT ^ 
dEoir) 



dT 



T 

ASoiT'), 



(12) 



where A5o(r') = AUo{T')/T' is the entropy jump 
obtained without any T-dependent parameters in the 
Hamiltonian, but where the quantity is to be evaluated 
at the temperature T' = T/Eq{t). Note that the pref- 
actor relating AS" to ASq would always be 1 irrespective 
of what Eq(t) is, if we had not included the contribution 
-T < dH/dT > to U. Using Eo{t) = 1 - t^, we find 



AS'(r) 



A5o(T'), 



(13) 



where H is an effective T-dependcnt Hamiltonian, 



precisely as in Ref.E3. Note the diffeipnce in the argu- 
ments of AS{T) and AS'o(T'). Ref.El concludes that 
within a line- liquid model with moderate values of ASq, 
substantially enhanced values for AS* are obtained, par- 
ticularly in the low-field regime, in agreement with ex- 
periments. The main factor in the enhancement is the 
denominator 1 — r^, which vanishes as T — > Tcmf- 

Note that the above procedure of substituting i7o with 
Eq{t) Hq does not in itself in any way assume that the 
physics of the vortex system in the low-field regime is 
determined exclusively by field.-induced flux lines. How- 
ever, were we to follow Ref.L3 and in addition assume 
that in the low-field regime, there only exists one rel- 
evant length scale in the problem, namely the magnetic 
length ao ~ 1/Vb, we wouldhe assuming that only field- 
induced vortices are relevant degrees of freedom on the 
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melting line. Our main point is that this may be ques- 
tionable in the low-field regime, and we will therefore 
refrain from utilizing such an assumption. 

If we insist on comparing-^iS'o(r') with results ob- 
tained using Ho, and not then we must fix T' to 
values obtained for the melting line in such calculations. 
Thus, T cannot vary arbitrarily between and 1, while 
fixing ASq independently. Rather, r and T' are related 
via T' = T/Eq{t). In calculations of A^o using Ho, we 
must therefore have T'/Tcmf < 1- Using Eo{t) = 1 — r^, 
we find T < (^5- l)/2. This gives enhancement factors 
(1 + r^)/(l - r^) < \/5 within the two-fluid model. If we 
express Eq. |l^ in terms of r' = T'/Tc, we obtain 

AS = Vl + 4t'2 ASoiT'); r' G [0, 1 > . (14) 

Similar enhancement factors may be found using the sim- 
plest mean- field ajpwroximation Eo{t) = 1 — r, as for in- 
stance used in Ref.E3. It would yield an enhancement fac- 
tor in Eq. |l^ given by 1+t'. In Fig. |ll| we have also plot- 
ted the entropy jump as obtained using a T-dependent 
prefactor in the Hamiltonian. We have used the results 
obtained using Hq and enhanced them by the prefactor 
in Eq. The inset of the figure shows the enhancement 
factor on the melting line obtained in our simulations. It 
varies quite slowly as a function of r' in the entire inter- 
val. Hence, even if we include the effect of Eo{t) on AS", 
we obtain an essentially field-independent entropy jump 
in the field regime considered in Fig. |ll|. For specificity, 
we have chosen Eo{t) = 1 — r^, and ignored the differ- 
ence between Tc and the mean-field ctitical temperature. 
We note also in this context that RefE3 finds an entropy 
jump of the magnitude we have found here within the 
lowest Landau-level approximation. Furthermore, Ref.E2l 
finds similar results using the isotropic Xy-model with 
/ = 1/6, iuragreement earlier simulations on the same fill- 
ing fractionlla. Note the large difference in fiUingiactions 
between the present work and the work of Refs.EBEa. For 
/ = 1/6, commensuration effects due to the numerical 
lattice are severe, and could conceivably lead to overesti- 
mates of the magnitude of AS'. This has been part of the 
motivation for pushing the simulations to the low filling 
fractions used in this paper. 

E. B-T Phase Diagram 

To estimate the real magnetic field induction B cor- 
responding to the flux-line densities considered in this 
paper, we use Eq.|| and take ^ab = 12—15 A. With this 
value of ^ab, we find the magnetic field corresponding 
to the smallest fiux-line densities considered (f=l/128) 
to be approximate 5 — 7T. In Fig. |l^ we show the 
f-T phase diagram originating from simulations of the 
XY model. The fiux-line densities / considered are 
1/f = 12, 14, 16, 20, 25, 32, 48, 64, 72, 84, 96, 112, 128. We 
see that the overall behavior of this phase diagram is con- 
sistent with the phase diagram in YBCO measured by 



SchiUing et alM, and Junod et alM. The FLL melting 
line at T^(f) separates the superconducting Abrikosov 
FLL phase from the incoherent vortex liquid phase, the 
latter being characterized by finite resistivity and strong 
diamagnetic fluctuations, with simultaneous loss of Bragg 
peaks in the FLL structure factor, flux-line integrity, and 
global phase-coherence in all directions. The remains of 
the zero field vortex loop "blow out" around Tbc2(/) de- 
stroys phase-coherence on all length scales, and thus sep- 
arates the incoherent vortex liquid phase from the normal 
metal phase. The melting line T„(/) decreases with de- 
creasing / with a positive curvature. Note also that the 
width of the critical region is large enough to influence 
the FLL melting transition over a sizeable field range. 
This field range is seen to extend up to / « 1/256 which 
we may conservatively estimate to be at least of order 
0- IT. 



V. CONCLUSION 

In this paper, we have investigated characteristics of 
the molten phase of the Abrikosov flux-line lattice via 
Monte-Carlo simulations on the three-dimensional uni- 
formly frustrated XY-mode\. Bragg-peaks in the static 
structure factor and phase-coherence along the direction 
of the applied magnetic field are both lost simultaneously, 
rendering the vortex liquid phase- incoherent. This be- 
havior is triggered by thermal excitations of closed vor- 
tex loops of diameters of the order of the average dis- 
tance between flux lines in the low-temperature lattice 
phase. On the melting line, this mechanism suffices to 
produce highly nontrivial vortex configurations with ap- 
preciable statistical weight on the template of field in- 
duced vortices. These configurations are characterized 
by a "percolation" of closed vortex-loops threading the 
entire sample in any direction. In particular, this is the 
case for directions transverse to the direction of the ap- 
plied magnetic field, which is tantamount to a loss of line 
tension of the field-induced flux lines. It renders a pic- 
ture of the molten phase of the flux-line lattice in terms 
of a liquid of well-defined, separated, and directed line- 
objects, invalid. Equivalently, a picture in terms of world- 
lines of 2D non-relativistic superfluid bosons is invalid in 
the liquid phase. An effective theory of the flux-line lat- 
tice melting and the vortex-liquid phase thus appears to 
present a formidable challenge involving the solution of 
a self-consistent coupled theory of field-induced flux-iine 
objects, and thermally induced closed vortex-loopsQEl. 
This coupling must evidently render the flux-line tension 
equal to zero in the liquid phase. Unfortunately, it is 
therefore doubtful that the intuitively appealing physics 
of directed polymers is particularly relevant for the vor- 
tex liquid phase. 

Scaling functions for the speciflc heat are calculated, both 
in zero and finite magnetic field. The zero-field results 
yield a sizeable critical region \T — Tc\/Tc « 0.1, corrob- 
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orating the notion that critical fluctuations of extreme 
type-II superconductors, i.e. vortex loops, will influence 
such phenomena as flux-line lattice melting over an ap- 
preciable range of magnetic inductions, possibly up to 
fields of order IT in moderately anisotropic supercon- 
ductors. The field range will depend on mass-anisotropy, 
since the width of the critical region and the low-field 
shape of the melting curve both appear to be influenced 
by the layeredness of the superconductor. 
The finite-field results for the scaling functions for the 
specific heat, as well as the obtained phase-diagram for 
an anisotropy parameter F = 3, are consistent with ex- 
periments on the slightly more anisotropic cuprate high- 
Tc superconductor YBCO, with F « 7. 
Finally, we note that columnar defects will not be partic- 
ularly efficient in enhancing the critical current density in 
a superconductor where the FLL melting line is strongly 
influenced by thermally excited closed vortex loops. (The 
influence of columnar defects on the vortex-sjstem was 
studied using Monte-Carlo simulations in RefO for a fill- 
ing fraction / = 1/2). The vortex- loop susceptibility 
should be sensitive to the phase-stiffness of the super- 
conductor. The phase-stiffness is in turn largely con- 
trolled by the superfluid density, and therefore also by 
the charge-carrier density. In order to avoid the detri- 
mental effects on transport properties in high- Tc super- 
conductors from a vortex-loop "blowout" , an increase of 
the charge-carrier density appears to be essential. 
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APPENDIX A: INTERNAL ENERGY 

In this appendix, we give a brief derivation of a gener- 
alized expression for the internal energy of a system with 
an effective T-dependent Hamiltonian. Consider a sys- 
tem in the canonical ensemble. For illustration, we will 
consider the well-known (P, V, T)-system. Our result for 
the internal energy U does not depend on the nature of 
the work-term. The system has a statistical distribution 
function given by the canonical law 



Z 



1 _H 

p——e "5", 
^ Z 



(Al) 



where the normalization constant Z is the canonical par- 
tition function 



E 

configurations 



(A2) 



where is a parameter of the distribution function which 
remains to be determined, such that 



configurations 



(A3) 



We insist that this normalization is to be maintained if 
the parameters V and are varied differentially. The 
Hamiltonian and hence the partition function will depend 
on V through the wall-potential of the problem. Let us, 
arbitrarily, write the partition function in the following 
way 

Z = e"*, (A4) 

where is a system-dependent parameter which also re- 
mains to be determined. When V V ^ dV and 
6 + d6, we will therefore also need to vary ^E* -I- d'^ 
in order to maintain correct normalization of p. Hence, 
we have 



E 



q>-~-H(V,6) 

e 5 



= 1 = 



■^ + d^-H(V+dV,9 + de) 
g 8+d9 _ 



(A5) 



Note that we have allowed H to depend on the statistical 
parameter 9. Expanding to first order in all differentials, 
we obtain 



E 



e ^ 



dH 



dO) 



dO 



1. 



(A6) 



Since the orginal distribution prior to changing V — > 
V -t- dV and 6 — > 6 + d9 also was normalized we ob- 
tain the following constraint on the differentials dV , dO, 
and 

Here, < .. > denotes a statistical average with respect to 
the original distribution function exp{{^ — H) / 9) . In or- 
der to make the connection to thermodynamics, we now 
compare the above with the "thermodynamic identity" 



dF dT^ 



PdV^—(F 
T ^ 



U)-P dV, (A8) 



where F ^ U — TS is Helmholz free energy, S is the 
entropy, and U is the internal energy. This comparison 
yields directly 



p ^ 


. dH 






dO 


dT ^ 




~T ^ 




F 


u = 


<H > 



= fcsT, 



(A9) 
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Note that 5* thus identified is the only choice consistent 
with F = -fesTlnZ. This then fixes U. Also, the ex- 
pression for U obtained in this fashion is identical to that 
obtained directly from the usual relation 



U 



dlnZ 
dp ' 



(AlO) 



with an assumed T-dependent Hamiltonian. 
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FIG. |l[ Specific lieat C per site, inplane structure factor S{\i±_ = 27r/5,7r/4), lielicity modulus along z-axis Tz, and helicity 
modulus along x-axis as functions of temperature for the system with vortex line density / = 1/20. The inplane structure 
function 5(kx) jumps discontinuously from 0.2 to precisely at Tm = 0.531Jx indicating that the FLL melts in a first order 
phase transition. At the same temperature, also shows a discontinuity from 0.6 to 0, indicating that the FLL melts directly 
into the incoherent vortex liquid with no global phase coherence along the applied magnetic field direction. At temperatures 
above Tm there is no global phase coherence in any direction. The specific heat also shows a (5-function anomaly precisely Tm. 
The broad specific heat anomaly at Tbc2 ~ 1.05Jx represents the remains of the zero-field Onsager vortex loop blowout. Note 
that for temperatures Tm <T < Tbc2 local superconducting phase coherence still exist, giving strong diamagnetic fluctuations 
in the liquid phase. The FLL depins from the numerical lattice at Td <^ Tm where vanish. Thus, the FLL melting transition 
at Tm S> Td is not affected by the numerical lattice. 

FIG. 1^. Intensity plots of the structure function S'(kx) for various temperatures for the system with flux-line density 
/ — 1/20. kx G [— 7r,7r] and ky G [— tt, tt] is along the horizontal and the vertical direction, respectively. The brightness in the 
plots is a measure of the magnitude of S'(kx). To enhance features we put all points where S'(kx) < 0.01 (noise level) to black 
and all points where ^(kx) > 0.05 to white. Precisely at Tm, the sharp Bragg-peaks in ^(kx) are converted into Bragg-rings, 
characteristic of a liquid. Thus, the FLL melts into a vortex liquid within a temperature region of AT — O.OOlJx- 

FIG. ^ Internal energy per site _E as a function of temperature for the system with vortex line density / = 1/20. The data 
are obtained from a cooling sequence using 3000000 sweeps per temperature. The internal energy has a discontinuous jump at 
Tm indicating a first order transition from an ordered state (FLL) to a disordered state (phase- incoherent vortex liquid) . This 
jump in the internal energy is used to determine the latent heat (entropy jump) at the FLL melting transition. The jump in 
E here corresponds to a jump in the entropy per vortex line per layer AS{f = 1/20) — OAks- 

FIG. ^. Snapshots of the vortex configuration for the system with vortex line density / — 1/20 for four temperatures, 
T/Jx = 0.26,0.50,0.54,0.70. For clarification we have shown only a part of the system; x,y € [0 : 20] and z € [0 : 40]. For 
T — 0.26Jx <S Tm, the fiux lines form a hexagonal lattice. Although there are many thermally induced defects attached to 
each flux line, they are nonetheless well defined quantities. For T = 0.50Jx ~ Tm, the FLL is still intact. Although the flux 
lines now contain many larger defects, they are still well defined. For T = 0.54Jx ~ Tm, the FLL has melted. For T ~ 0.54Jx, 
it is seen that the flux lines are no longer well defined quantities. There exists at least one way for a flux -line to thread the 
system in any direction. For any vortex-conflguration, therefore, there exists at least one flux line threading the sample in the 
direction perpendicular to the magnetic field. 

FIG. 1^. The inplane structure function ^(kx) as a function of temperature for several vortex line densities /. For a given 
/, 5(kx) with the corresponding value of kx shows a sharp drop from ~ 0.2 to at a well-defined FLL metling temperature 

Tmif). 

FIG. 1^. The helicity modulus along the field direction as a function of temperature for several flux-line densities /. For 
all densities, shows a sharp drop towards zero precisely at the corresponding FLL melting temperature Tm{f)- Thus, the 
FLL melts directly into an incoherent vortex liquid. 

FIG. 1^. Monte-Carlo results for the specific heat per site of the anisotropic 3-DXy-model as a function of temperature for 
several vortex line densities /. The system sizes depends on filling fraction, as explained in the text, and F = 3. For clarity the 
nth curves is shifted by an amount 0.2 * n upwards. For each / there is a spike at a /-dependent critical temperature Tm{f) 
indicating a flrst order phase transition. 
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FIG. g. Monte-Carlo results for the specific heat of the anisotropic 3DXY-mode\ as a function of temperature for several 
flux-line densities. For decreasing / (decreasing magnetic field induction B) the crossover temperature TBc2{f) increases and 
moves towards the zero field critical temperature TBc2{,f ~ 0) ~ Tc- The broad anomaly (cusp) in the specific heat becomes 
sharper and the maximum height of the cusp increases. Thus, for decreasing / the specific heat evolves smoothly to the zero 
field specific heat singularity at Tc- The spike in the specific heat at Tm{f) <C TBc2{f) each graph is hidden in the noise of the 
other graphs and is therefore hard to recognize in this particular figure. While in zero magnetic field the vortex loop blowout 
is the mechanism for the second order phase transition at Tc, in finite magnetic field the vortex loop blowout at TBc2{f) is only 
a crossover. The phase transition in systems with finite vortex line densities take place at a lower temperature, Tm(/), where 
the vortex line lattice melts. 



FIG. |9|. Monte-Carlo results for the specific heat of the anisotropic 3-DXy-model in zero magnetic field, for various system 
sizes L X L X L with L = 32, 48, 64, 72, 96, and two values of the anisotropy, F = 1 and F = 3, scaled according to Eq. ^. Here, 
t — {T — Tc)/Tc. The region of data collapse gives the width of the critical region. Note that this region is slightly for F = 3 
than for F = 1. 



FIG. Monte-Carlo results for the specific heat of the anisotropic 3DXY-mode\ for a number of filling fractions / given 
by 1// = 12, .., 128 with corresponding system sizes as explained in text, and anisotropy F = 3, scaled according to Eq. ^ The 
results are in good agreement with the experimental results of Schilling et al, and Junod et al. 

FIG. |ll| The entropy jump per vortex line per layer A5'(/) at the FLL melting transition for several vortex line densities 
/. The filled circles represent the results obtained with a T-independent Hamiltonian, Eq. |l| The open circles represent the 
results obtained including a T-dependent prefactor in the Hamiltonian. We see that AS{f) essentially does not depend on / 
in this regime of filling fractions /, regardless of whether T-dependent prefactors are included in the Hamiltonian or not. The 
inset shows the enhancement factor in Eq. |l^. 

FIG. |l^. The helicity modulus perpendicular to the field direction as a function of temperature for several vortex 
line densities /. For each /, Tx vanishes at a temperature Td{f) significantly lower than the corresponding FLL melting 
temperature Tm{f)- The artificial pinning potential of the numerical lattice therefore does not affect the FLL melting transition 
at T^if) » Taif). 



FIG. g_3[ The f-T phase diagram for the uniformly frustrated 3d XY model. The applied field is along the crystal c-axis, 
the anisotropy parameter F = 3. The FLL exhibits global phase coherence along the applied field direction. The FLL phase 
is separated from the incoherent fiux-line liquid phase by the melting line Tm(/). The melting transition is a first order phase 
transition with an entropy jump AS{f) ~ O.lks for the anisotropy F = 3 and field regime considered in this paper. In the 
incoherent vortex liquid phase Tm{f) < T < TBc2{f), there is only local, but no global, phase coherence in any direction. 
At finite fields, between the incoherent vortex liquid and the normal metal phase, there exists a broad crossover region where 
a blowout of thermally induced closed vortex loops takes place, eventually also destroying superconductivity on short length 
scales. The width of the crossover regime is obtained from scaling behavior of the specific heat. Another, consistent, method 
of obtaining this width, is to estimate the temperature regime which correponds to an uncertainty of 10% in the maximum 
value of the specific heat anomaly at TBc2{f)- Since this anomaly becomes broader with increasing field, the cross-over region 
becomes wider. This is also confirmed from the scaling results for the specific heat. 
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